Study of the Distillability of Werner States Using Entanglement 
Witnesses and Robust Semidefinite Programs 

Reinaldo 0. Vianna^* and Andrew C. Doherty^ 
*t School of Physical Sciences, The University of Queensland, Queensland 4 072, Australia 
* Departamento de Fisica - ICEX - Universidade Federal de Minas Gerais 
Av. Antonio Carlos 6627 - Belo Horizonte - MG - Brazil - 31270-901 

(Dated: February 1, 2008) 

Abstract 

We use Robust Semidefinite Programs and Entanglement Witnesses to study the distillability of 
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of the state space which was previously conjectured to be undistillable. We also introduce bases 
which yield interesting expressions for the distillability witnesses and for a tensor product of Werner 
states with arbitrary number of copies. 
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I. INTRODUCTION 



Maximally entangled states are the main resource in Quantum Information (QI) process- 
ing. Protocols like teleportation of quantum states [1] and entanglement based quantum 
cryptography [2], just to cite the two most emblematic, work with perfect fidelity only when 
maximally entangled bipartite states are available. ( As maximally entangled bipartite states 
are equivalent by local unitary transformations, they are usually referred to as singlets in the 
jargon of QI.) On the other hand, even if one has a source of perfect singlets, the ever present 
decoherence, due to interactions with the environment, degrades these states to mixed form 
with reduced entanglement. Bennett et al. [3] showed that this practical difficulty could 
be circumvented by distilling singlets from mixed states. This process involves only Local 
Quantum Operations and Classical Communication (LOCC) and, in principle, is able to 
purify any mixed state to a singlet form, given that an arbitrary supply of the former is 
available. The process works at the expense of sacrificing many of the mixed states in order 
to concentrate their entanglement to a singlet. Horodecki et al. showed [4,5] that only 
states which violate the Peres criterion [6] can be distilled, i.e., the non-positivity of the 
partial transpose (NPT) is a necessary condition for distillability. It was then realized that 
there are entangled states which are not directly useful in QI processing, they are said to be 
bound entangled [4], and the states with positive partial transpose (PPT) are of this kind. 
Nevertheless, these states can be activated [7], in the sense that, used in conjunction with 
NPT states, they can enhance the fidelity of teleportation. Therefore we have two kinds of 
entanglement in nature, namely, bound and free. The set of bound entangled states includes 
all the PPT ones, but it is not known if there are NPT states in this set. It was conjectured 
by DiVincenzo et al. [8] and Dur et al. [9] that, in fact, there exist bound entangled NPT 
states. 

It can be shown that any bipartite NPT state can be transformed by LOCC to a Werner 
state [10], keeping the fidelity to the singlet. The process is performed by twirling (see [9], 
for example) the states through the action of bi-local unitary operations (U <g) U). Thus the 
study of distillability of arbitrary bipartite states is reduced to the distillability properties 
of Werner states. 

Formally, a bipartite state (p G B(7i^ <8> H-b)) is distillable if and only if there exists a 
Schmidt rank two pure state (|^) G Hj, <8> H-b), in the Hilbert space in which p acts, such 
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that (\l/|(/y g,Ar ) TA |\I/) is less than zero for some finite integer N [5,8,9], Ta meaning partial 
transposition. When this condition is verified for some N, p is said to be iV-distillable. In 
particular, all the bipartite entangled states of the kind qubit-qudit (2 x N) are 1-distillable 
[5,11]. 

In the same fashion that the entanglement of a state can be decided by an Entanglement 
Witness operator [11], Kraus et al. showed that the distillability also can be decided by 
means of a kind of witness operator [12]. In this paper, we show how to calculate these 
distillability witnesses using Robust Semidefinite Programs (RSDP) [13] and apply it to 
study one- and two-distillability of Werner states. Starting with some definitions in section 
II, we revise the RSDP formalism in section III. Sections IV and V present numerical results 
for the distillability of Werner states in the one- and two-copy cases. In section VI, we derive 
some interesting expressions for the distillability witness and for a tensor product of Werner 
states, with arbitrary N, and then we conclude. 

II. DEFINITIONS 

The set of non-entangled (separable) states is convex and closed, therefore it follows from 
the Separating Hyperplane Theorem that there exists a linear functional ( hyperplane) that 
separates an entangled state from this set, this results in an Entanglement Witness [11]. 
Thus an EW (W) is a Hermitian operator with non-negative expectation value for all the 
separable states, but which can have a negative expectation value for an entangled state, in 
this case, the state is said to be detected by the EW. An EW which can be written in the 
form 

W = P + Q T 1 A + Q T 2 B + ... + Q T N z 1 (1) 

with P and Qi positive operators, is said to be decomposable, and it is non-decomposable if 
it cannot be put in this form. Only non-decomposable EWs can detect PPT states. When 
the EW (hyperplane) is tangent to the separable set, it is said to be optimal (see [14], for 
example). 

To be distillable [4,12], a bipartite state p G B(Ha <8> Hb), or a finite tensor product of it 
(p 0jv ), must have a projection which is NPT on a four-dimensional subspace, that is, given 
a Schmidt rank two state, 

|*) = si|ei/i)+s 2 |e 2 / 2 ), (2) 
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where ei,e2 and /i,/2 are bases defining bi-dimensional subspaces in H.j^ N and 1-Lb® n , 
respectively, and si, S2 are the Schmidt coefficients; p is distillable if the inequality 

(q\(p® N ) TA \q) < (3) 

is satisfied for some arbitrary |\&), of the form of (2) , and some finite integer N. 
We use the following parametrization of the Werner states [10] : 

pw ~ d^ + dp' [ } 

with — 1 < (3 < 1. p w is separable for (3 > — ^ and 1-distillable for (3 < — \. F d is a swap 
operator for two qudits, 

d 

^ = Z>'>o«i, (s) 

and its partial transpose is the bipartite maximally entangled state: 

^ = ^ a = ^X>>0jI- ( 6 ) 

Id is the identity in the space of the two qudits (Tr(I d ) = d 2 ). 

Equivalently to the inequality (3), Kraus et.al. showed [12] that the distillability of an 
arbitrary state p (we will consider only the bipartite case p G BiTi^ <g> H-b)) can be decided 
through the operator: 

W N = P 2 ® (p T -f N , (7) 

with Wn acting in (H 2A <g> 7Y 2 b) <8> (H-a ® 7~Lb)® N , 'H2A (^2b) being the Hilbert space of a 
qubit belonging to A (B). If Wn is not an EW, then p® N is N-distillable. If Wn is a non- 
decomposable EW, then the PPT entangled state it detects (tt) activates p 0Ar , i.e., p 0Ar ®7r 
is 1-distillable. When W 7 ^ is decomposable (and in the case of Werner states it happens to 
be a positive operator), then p® N is undistillable and unactivable. When Wn happens to be 
an EW, Wl A and W T N B are optimal EWs. 

III. OPTIMAL WITNESSES VIA ROBUST SEMIDEFINITE PROGRAMS 

Given a bipartite state p G B(1~La <E> Kb), we want to determine its optimal entanglement 
witness W p . We will use the method introduced in [14], which we briefly describe in the 
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sequel. Let E be the set of separable states in 13 (Ha <8> H B ). We define the following set 
(W) of entanglement witnesses: 

W = {W / W e B(H A ®H B );W ] = W;Tr(Wa) > 0,Va G ^;Tr(W) = 1}. (8) 

W p is defined by: 

mm Tr(Wp)=Tr(W p p), (9) 
and can be determined through the following RSDP: 



min Tr(Wp), 
w 

f 

Wt = W 

subject to { Tr(W) = 1 (10) 
(^a\W\^ a ) > V|Va) e ?U 
This is a NP-HARD problem and, in practice, we relax it to a Semidefinite Program (SDP) 
by taking a finite number (n) of pure states (\iPa)) to represent the whole Hilbert space 
H.A- Thus we replace an infinite number of constraints ( cf. last line of (10)) by a finite 
set. If dim(H a) is d, the kets \ip l A ) can be chosen as an uniformly distributed sample of unit 
complex vectors (q) in C d , and with infinite n this program would yield the exact witness 
W p . In our calculations, we use an interior point algorithm to solve the SDP [15]. 

In [16], it was shown that W p yields the random robustness (Rr = —Tr(I)Tr(W p p)) 
of p, i.e., the minimal amount of mixing with the identity necessary to wash out all the 
entanglement. Thus, the state a = (p + Rrl /Tr(I))/(l + Rr) is in the border between 
separable and entangled states. 

Our main goal is to decide if the operator W N (cf. (7)) is an EW. We will do so by 
determining a state for which Wn could be an optimal EW, in the sense of (9). If we find 
such a state, we compare its optimal EW with the expression of Wn, and this will tell if 
Wn is or is not an EW. 

IV. ONE-COPY CASE 



We will apply the RSDP techniques to calculate optimal EWs (cf. (10) to investigate 
the distillability properties of Werner states in the one-copy case. We will show that the 
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distillability is related to the properties of an EW for a certain family of PPT states. 
We want to know if the Hermitian operator Wi(f3) G B(H_a ® Hb) defined by: 

W 1 (0) = P2®f&, (11) 

is an entanglement witness. 

We will show that, for — | < f3 < — |, this operator is indeed a witness and, for (3 = 
— 1, it is the optimal witness (W w ) for a certain family of PPT entangled states (ir). Our 
numerical calculations will be restricted to qutrit Werner states(d = 3), therefore dim{Ti.js) = 
dim(H B ) = 2xd = 6. 

We introduce the orthogonal basis: 

B l = P 2 ® P 3 ; 
B 2 = P 2 ®(h-P,)- 
B 3 = (I 2 - P 2 ) (g) P 3 ; 
, B 4 = (I 2 -P 2 )®(I 3 -P 3 ). 

Note that B i /Tr(B i ) is an entangled state in B(Ha <S> T~Cb), an d only B 4 is PPT. In particular, 
B 1 is the maximally entangled state. 
Wi(f3) can be recast as: 

W M = d^Vpd^ 1 + pd)Bl + B2] - (13) 

In this basis, the states can be written as: 

P = TLiPiBi/Tr(B i ), 

T,t=iPi = i' Vi> o. 

Note that the state space is a 3-dimensional polytope defined by the constraint J2t=i Pi = 
1, pi > 0. Optimal witnesses for these states have the form: 

W = Zt 1 c i B t /Tr(B t ), 

Ej=l °i = !> C i G K - 

Assuming Wi(/3) is a witness for — 1 < (3 < — |, we ask for the PPT state (7r) it detects the 
most. It is done through the following SDP: 



min7V[Wi(/3)7r], 
6 



subject to < 



7T > 
Tr(vr) = 1 

7T Ta > 



(16) 



We observe that the optimal solution (re*) of the above SDP is independent of (3 and mini- 
mizes this other SDP: 



minp 

' 7T = (1 - p)B 1 /Tr(B 1 ) + pB 4 /Tr(B 4 ) 

= 7T 

subject to \ 7r > (17) 
Tr(vr) = 1 
tt Ta > 

The optimal p is 0.8571, yielding the state ir*. The optimal EW for this state, obtained by 
means of the RSDP (10), furnishes Tr{W n .ir*) = -0.0095. 
Comparing Tr[Wi((3)%*] with Tr(W n *ir*), we observe that 

TriW^ir*)] < Tr(Ww*), V - 1 < (3 < - X -. (18) 

These calculations are sufficient to show that W\{f3) is not a witness for (3 < — |, for 
it gives an expectation value that is lower than that of the optimal EW. On the other 
hand, for -\ < (3 < -|, we observe that 7V[Wi(/3)7r*)] > Tr(W n *ir*), and for f3 = -|, 
Tr[Wi( — |)7r*)] = Tr(W n *7r*). With W w * written in the form (15), our calculations converge 
to ci = — -tt, C2 = ^| and C3 = C4 = 0, which are the parameters of W\{— \). Therefore 
Wi{— \) is the optimal witness for 7r*. This result can be obtained using a large sample of 
random \tp l A ) ( cf. (10)) or through the following deterministic recipe. 
Consider the state (a) defined by: 



7r* — Tr[Wi( — ^)ir*]I 



(19) 



1 - (2d) 2 Tr[W 1 (-^)n*] 

If Wn* and W\(— \) coincide, they yield the random robustness of it* [16], and a is a state 
in the border between separable and entangled states (viz. Sec. III). Therefore a contains 



information about the border of the separable set. For each eigenvector (l^fc)) of a, we form 
the state p\ = Tr B (\^k)(^k\)- Then we use the eigenvectors of the {p k A } as a sample of 
states for the SDP (10). In this case, this recipe yields 216 states, but just 24 already yield 
the exact result. 

Now we want to show that W\{($) is a witness for — | < (5 < — |. For (3 = — |, we know it 
is an optimal witness. Notice that W\{— |) is a positive operator. Looking up (13), it is easy 
to see that, for any state a and for p 1 < f3 2 , Tr\Wi(Pi)<j)] < Tr[Wi(f3 2 )a]. In particular, 
7V[Wi(-D<7] < Tr[Wi(/3 2 )o-]. If a is a separable state, 7V[Wi(-±)<7] > and therefore 
Tr[W±(P > —\)<j] > 0, showing it is an entanglement witness. 

All the calculations we have done can be understood more easily by means of Figs.l and 
2. Fig.l shows a 2-dimensional projection of a 3-dimensional plot of the state space. This 
picture is obtained as follows. We randomly draw 10 6 states p. Out of each p, we build a 
border separable state (a), and a state (0) in the hyperplane defined by Wi(f3), namely: 

= p-Tr(W pP )I . 

l-36Tr(W» ' (20) 

V l-36Tr[Wi(/3)p] " 

These states are rewritten in the zero trace basis (J, G\, G 2 , G 3 ), 

G\ = 8Bx — B 2 ; 

G 2 = 8B 3 - B 4 ; (21) 
G 3 = -3(Bi + B 2 ) + B 3 + B 4 ; 

and their coefficients are plotted. We clearly see that the planes Wi((3) have a common axis, 
which is parallel to G3. In the picture, we can also see the state 71* ( cf. (16, 17)), which 
is in the plane W\{— |). Notice that the plane Wi(— |) is tangent to a face of the polytope 
defined by the separable states. It can be clearly seen in Fig.2, which is the 3-dimensional 
picture. Fig.2 also clearly illustrates the concept of Optimal Entanglement Witness, i.e., 
a hyperplane tangent to the separable set. Notice that in Fig.l, we have a family of non- 
optimal EWs, the planes for — | < f3 < — |. Therefore, these calculations show that the 
Werner states are 1-undistillable for — | < /3 < — |. 

Of course the 1-undistillability of Werner states is not a novelty [8,9]. The interesting 
result here is the technique to decide if the Kraus-Lewenstein-Cirac operator (7) [12] is an 
EW. The strategy was to compare the candidate to EW with the optimal EW of a certain 
PPT state, which can be obtained with arbitrary precision by means of the RSDP (10), and 
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FIG. 1: A projection of the state space showing the set of separable states (just border states - 
cf. (20)), the planes W\{(3) and the PPT state ir for which W\{— |) is the optimal witness. The 
planes separate the state tt from the separable set. 

to show that the candidate operator converges to this EW. In this sense our calculations 
are exact, leaving no room to doubt if the candidate operator is or is not an EW. This 
technique extends straightforwardly to higher dimensions. The other interesting result is 
that we obtain a family of PPT entangled states, vr(p), 0.8571 < p < 1 (cf. (17)), which 
activate the Werner states in the interval — \ < [3 < — ~, i.e., p w <g> n(p) is 1-distillable. In 
[17], similar results about the activation of Werner states were also obtained. 
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G2 



-2 -2 



G1 



FIG. 2: A three dimensional picture of the state space showing the set of separable states (just 
border states - cf. (20)), which is a polytope, and the plane W\{— |) sitting on the polytope. W\{— ^) 
is an Optimal Entanglement Witness, i.e., a plane which is tangent to the separable set. 

V. TWO-COPY CASE 

Now we will apply the techniques we have developed in the one-copy case to study the 
distillability of Werner states in the two-copy case. We will determine the optimal EW, by 
means of the RSDP (10), for a family of PPT states. This will show that the Werner states 
which are 1-undistillable are also 2-undistillable. 

The calculations for the two-copy case mirror the one-copy case and we arrive at analogous 
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conclusions. The orthogonal basis to expand the witnesses and states is: 



B 1 


= P2 


g>P 3 <S 


>P*; 






B 2 


= P2 


8>[(/ 3 


-P 3 ) 


g> p 3 + p 3 ® (J 3 - 


^3)]; 


B 3 


= P2 


2>(h- 




Hh-Ps)}; 




B, 


= (h 


-P2) 


®p 3 6 


5P 3 ; 




B 5 


= {h 


-P2) 


®[(/ 3 


- P 3 ) <g> P 3 + P? e 


Hh-Ps)}; 


B e 


= {h 


-P2) 




- p 3 ) ® (/ 3 - p 3 )] 





(22) 



Again, if normalized, these are entangled states, with B 1 the maximally entangled state, 
and B 6 the only PPT state. The state space is a 5-dimensional polytope defined by: 



W 2 {(3) = P 2 ® {p T w A )® 2 reads: 
W 2 {(3) = (- 



P = HiftS i /Tr(S i ); 
Ei=iP» = ! ? Pi > 0- 



-) 2 [(1 + /?d) 2 Bi + (1 + (3d)B 2 + S 3 ]. 



(23) 



Using a SDP analogous to (16), we obtain that the PPT state (tt) most detected by W 2 ((3) 
has the coefficients (0.0278, 0.2222, 0, 0, 0.0833, 0.6667). We also obtain that W 2 {-\) = W n 
and 7V(W» = -0.0019. We can derive different families of PPT states for which W 2 {-\) 
is optimal. A particularly interesting one reads: 

B 2 Pf; 



7T 



+ p- 



(25) 



Tr(P 2 ) "Tr(P 6 )' 

with p exactly the same as in the one-copy case, namely, 0.8571. For this state we have 
Tr[W 2 (-±)ir} = -0.0013. 

As our calculations show that W 2 (— \) is an optimal witness, W 2 (f3) is an entanglement 
witness for — \ < (3 < — |. Notice that these calculations are exact and show that qutrit 
Werner states are 2-undistillable in this interval of (3. On the other hand, they can be 
activated by the the families of PPT entangled states detected by W 2 (— |). The best results 
so far showed 2-undistillability analytically [8,9] in the region -0.417 < 13 < -\, and 
provided numerical evidence in — \ < (3 < —0.417. 

We note that our calculations have been made in a laptop with 1 GBytes of memory, and 
we used the symmetry of the Werner states to reduce the parameters in the optimization 
problem. In a larger computer, maybe the 3-copy case could be handled, but the 4-copy case 
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would need more than 5 GBytes, just to load the basis set (28). One could try to explore 
the symmetry to reduce the matrices size, but in face of the constraints in the Robust SDP 
( cf. last line of 10), which is also the most memory consuming part of the calculations, it 
is far from trivial. 



VI. THE W N OPERATOR 

Although we are computationally limited to calculations for the 2-copy case, we would like 
to understand the properties of the operator when the number of copies (A) increases, 
hoping to shed light on the general problem. Indexing the copies, we can write explicitly: 

W N = P 2 ®pli®pl A 2 ® • • -®pZ a n . (26) 

Note that it is normalized (Tr(W N ) = 1) and tracing out the Ath copy yields the operator 
for A — 1 copies: 

Tr N (W N ) = W N -l (27) 

Now we introduce a basis set which allows to write Wn as a polynomial, generalizing equa- 
tions (13) and (24). 

Define the following basis of orthogonal projectors for the A-copy case: 

B? = P 2 ®Pf N - 

Bj+i = P2® ^ np! [N ~ j) ® (/„ - p d n (28) 

B% +1 = P 2 ® (h - P d r N ; 

with Bf +1 G 13 (Ha <8> Hb)- The Pi permute the elements in the tensor product, yielding an 
expression which is totally symmetric under exchange of any P d and (I d — P d ). In this basis, 
Wn has the diagonal representation: 

W N = J2f=o^3+i B j+n 



A 3+ 1 ~ (cP+ d p) N ■ 

We will show the correctness of this expression by induction. Note first that 

,N v fN\ /J2 



(29) 



Tr(B» +1 )=^)(d'-iy, (30) 

where (^) is the binomial coefficient. The basis for A + 1 copies is related to the A-copy 
basis through the recurrence relation: 

Bf +l ®P d + Bf ® (I d - P d ) = B?+\ (31) 
12 



If we normalize the basis (28): 



bi+1 -0(<p-iy (32J 



we can rewrite (29) as: 

(33) 



W N 



Now it is easy to see that the trace of (33) is 1: 



(d 2 + dpy 



Tr{W N ) = ^ 3=0Ki ') M , ' V = I- (34) 



To finish the proof of the correctness of (29), we build W N+1 by appending p w to W N : 

w N+1 = w N ® P l A = 

WTW^ S^ 1 + W) N - S B? +1 ® [(I d - P d ) + (1 + d(3)P d \. 

Splitting this sum in two parts, and redefining the index in the second sum as j + 1 = k, we 
obtain: 



W N+l = (d , +d \ )N+i {Ef =0 (l + d/3) N+1 -*B? +1 ® P d 
+ EjLYU + d0) N+1 - k B» ® (I d - P d )}. 
Writing out explicitly the terms for j = and k — N + 1, we arrive at: 



(36) 



+ Ef=i(l + dp) N + x -*[B? +x ®P d + Bf ® (I d - P d )] (37) 
+£# +1 ®(/ d -P i )}. 

Finally, using the recurrence relation (31) and the basis definition (28), we obtain the desired 
result: 



W N+1 = jp+kpmW + dfl N+1 B™ + Ef =1 (l + d0) N+1 -*B>!ff + i?^ 1 } = 

Once the correctness of (29) is proved, we highlight an interesting property of the Wn 
operator for Werner states. If < \d 2 +d{3\, then all the eigenvalues of W N (f3) go to zero 

when N tends to infinity, for any — 1 < f3 < 1. This is an expected property. If Wn{— \) is an 
EW, for some N, it is necessarily optimal, and if properly normalized, it furnishes the random 
robustness (Rr) for a family of entangled states (ir), i.e., Rr = — {2d) 2N Tr{W^{—^)^) [16], 
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and we see that the entanglement, as measured by the random robustness, increases with 
N. 

We can also work out an expression analogous to (29) for a tensor product of Werner 
states: 



^ N _ (i d + (3F d r N 

{d 2 + dpy 

We construct the following basis set, which has the same algebraic structure of (28), 
although it is not orthogonal and only the last element is a positive operator: 

AN _ f®N. 
A l — Id ' 

A? +l = S[ff N -> ) ®(I d -f d r] ] (40) 
A N+1 = (I d - f d r N ; 

with A^ +1 e B(H A ® H B ), and f d = F d /d. Note that this basis is obtained by discarding P 2 
in (28) and by taking the partial transpose of P d . In this basis, p® N reads: 

1 N 

(«) 

Note that A% +1 is a fully separable positive operator. In particular, if we take (3 = — \ 
then 

AN AN 
n ®N _ S±N+1 _ S±N+1 (A0 \ 



Tr(A») {d*-iy 



Normalizing the basis: 
p m reads: 



AN AN 
N _ ^'+1 _ ^j+l 



N _ J + i _ J + l (AO\ 



Then, for /3 = — ^ , p 0Ar = ajy+i- When we take iV to infinity, the binomial coefficients in 
(44) with j 7^ and j ^ N dominate the sum, but nothing special seems to occur. 



VII. CONCLUSIONS 

We have done exact numerical calculations which show the 2-undistillability of qutrit 
Werner states in the region —\<(3<—\. We have shown that Wn{— \) is an optimal 
entanglement witness (for N = 1,2), and this is our certificate of 1- and 2-undistillability. 
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We have derived families of PPT entangled states which activate the 1- and 2-undistillable 
Werner states. We have introduced a basis of orthogonal projectors to expand Wn{@), which 
shows that this operator is a polynomial in (1 + d(3), and it acts in a state space which is a 
polytope. The eigenvalues of W^{f3) tend to zero, when we consider infinite many copies of 
Werner states and it is a property related to the random robustness of the states it detects. 
We have also introduced a basis which provides a simple polynomial expression for a tensor 
product of Werner states. 
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